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Abstract 


This  note  derives  the  distribution  of  electrical  spreading  currents 
along  the  length  of  solid  conducting  objects  for  which  the  length  sub¬ 
stantially  exceeds  the  width.  Sources  and  sinks  of  DC  (or  very  low 
frequency  AC)  current  are  placed  at  one  end  of  the  object,  and  the 
fall-off  of  spreading  current  is  studied  as  a  function  of  length.  The 
fall-off  can  be  great;  for  instance  in  the  case  of  a  solid  rectangular 
object,  the  fall-off  of  spreading  current  along  the  length  is  27dB  per 
unit  width.  Comparison  is  made  to  inductive  coupling,  which  becomes 
important  as  frequency  increases. 
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1  INTRODUCTION 


When  we  inject  a  current  into  a  solid  conducting  body  at  one  point 
on  its  surface,  and  allow  the  current  to  flow  back  out  at  another  point,  the 
current  distribution  within  the  body  spreads  out  around  the  direct  path 
from  source  to  sink.  This  is  called  “current  spreading” ,  and  is  important  in 
semiconductor  design,  where  crosstalk  among  different  devices  on  the  same 
chip  must  be  minimized  [1]. 

In  this  note  we  are  particularly  interested  in  the  decay  of  DC  current 
spreading  along  the  length  of  a  long  conducting  body  when  both  source 
and  sink  are  at  one  end.  We  will  find  that  the  magnitude  of  spreading 
current  generally  falls  off  like  exp{—kz),  where  2  is  the  distance  away  from 
the  current  source  along  the  long  axis  of  the  body,  and  k  is  a  characteristic 
decay  coefficient. 

In  the  language  of  electrical  engineering,  we  drive  the  body  as  a  (poor!) 
transmission  line  with  DC  in  differential  mode  at  one  end,  and  measure  the 
decay  along  the  distance  2  as  some  number  of  dB  per  unit  length, 

K  =  20\ogiQ{e)k  dB/length.  (1-1) 

Note  that  k  is  independent  of  the  resistivity  of  the  material  (as  long  as  the 
body  is  uniform  in  composition);  it  depends  only  on  geometry  of  the  body. 
The  current  source  is  assumed  to  supply  the  necessary  voltage  to  drive  the 
current;  the  decay  is  due  to  geometrical  effects  on  the  current  distribution, 
as  it  wends  its  way  through  the  resistive  material.  Experimentally,  we  can 
measure  the  fall-off  of  voltage  1^(2)  over  a  cross-section  of  the  body,  as  a 
function  of  distance  2  in  the  long  direction  of  the  body;  see  Figure  1.  In  fact 
this  measurement  would  be  easy  to  do  and  is  strongly  recommended  for  any 
practical  applications. 
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Figure  1;  Measurement  of  current  spreading  for  a  long  object.  We  inject  a 
known  current  in  differential  mode  at  the  end  of  the  semi-infinite  object.  We 
then  measure  voltages  V{z)  acro.ss  the  object  at  various  distances  2  in  the 
long  dimension.  The  body  has  cross-section  S  bounded  by  curve  C. 

At  DC,  the  only  important  length  scale  is  the  cro.s.s-.section  .size.  For  AC 
at  all  but  the  lowest  frequencies,  another  length  scale  becomes  important, 
namely  the  skin  depth;  and  inductive  effects  will  become  important  too.  We 
begin  by  .studying  the  DC  problem.  Finite-frequency  effects  including  .skin 
depth  and  induction  will  be  briefly  di.scu.s.sed  in  Section  4. 
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2  DC  EQUATIONS  AND  THEIR  SOLUTION 


We  use  cartesian  coordinates  (x,  y,  z)  and  model  our  solid  body  as  a 
right  cylinder  with  long  axis  along  2.  The  cross-section  of  the  body  is  a 
region  S  in  the  {x,  y)  plane  bounded  by  a  curve  C,  and  is  constant  in  2. 
One  end  of  the  body  lies  in  the  plane  2  =  0,  and  the  body  extends  over 
0  <  2  <  00,  and  so  is  half-infinite. 

Hollow  cylindrical  bodies  can  be  equally  w'ell  modelled  by  using  a  curve 
C  which  consists  of  two  (or  more)  separated  components,  one  inside  the 
other.  A  pair  of  separate  conducting  bodies  could  be  modelled  by  a  cross- 
section  S  consisting  of  two  separate  disconnected  components;  however  we 
will  take  S  to  be  connected  unless  otherwise  specified. 

The  body  is  composed  of  material  of  some  constant  resistivity  a.  We 
study  the  stationary  Maxwell  equations  [2]  in  which  all  time  derivatives  are 
neglected,  and  so  all  fields,  and  the  current,  will  be  functions  of  (x,  y,  2)  but 
not  time  Steady  electrical  current  flow  within  the  body  is  described  by 
the  current  density  J  and  obeys  the  equations 

V  •  J  =  0,  J  =  (tE  (2-1) 

where  E  is  the  electric  field.  On  the  sides  of  the  cylinder  no  current  flows 
across  the  surface,  so  the  boundary  condition  is 

n  •  J  =  0  on  C  for  all  2  >  0  (2-2) 

where  n  is  the  unit  normal  to  C  in  any  plane  2  =  const.  On  the  end  2  =  0 
of  the  cylinder  we  impose  some  source  Jo  of  injected  and  sunk  currents, 

z  •  J  =  Jo(x,  y)  at  2  =  0.  (2-3) 

Hn  this  section,  frequency  is  DC,  or  so  low  that  inductive  effects  can  be  neglected;  for 
the  latter  see  Section  4. 
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The  electric  field  can  be  expressed  in  term  of  the  gradient  of  the  electostatic 
potential 

E  =  -V$  (2-4) 

and  eliminating  both  E  and  J  from  Eqs.  (3. 2, 3. 5)  we  obtain  equations  deter¬ 
mining  $  a^ 


= 

0  inside  the  body 

(2-5) 

n  •  V<1>  = 

0  on  C  for  all  2  >  0 

(2-6) 

—a~^S{x,y)  at  2  =  0 

(2-7) 

dz 

which  describe  a  well  posed  boundary  value  problem  with  Neumann  bound¬ 
ary  conditions.  (For  notation  and  methods  see,  e.g.,  [2]). 


This  problem  can  be  solved  by  separation  of  variables, 


(I)(x,  y,  z) 

(2-8) 

=  0 

(2-9) 

n  •  Vip 

=  0  on  C 

(2-10) 

dz'^ 

(2-11) 

with  corresponding  boundary  conditions  at  2  =  0,  see  below,  where  A'  >  0  is 
a  separation  constant,  and  where  the  Laplacian  operator  in  the  (x,  y)  plane 
is 


V?  = 


5" 


+ 


(2-12) 


*  dx^  ' 

We  can  form  any  solution  $  by  use  of  the  eigenfunctions  of  the  transverse 
Neumann  boundary  value  problem, 


V'^t'il^iix.y)  = -\irpi{x,y),  n-V'ip  =  0  on  C 
with  eigenvalues  Aj,  i  =  0, 1, 2  . . .  by  means  of  the  expansion 

00 

^{x,y,z)  =  '^Cii;i{x,y)exp{-kiz) 

1=0 


(2-13) 


(2-14) 
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where  Aj  =  kf.  The  eigenvalues  Aj  will  be  assumed  ordered  by  nondecreasiug 
value,  and  the  eigenfunctions  'ip{x,  y)  will  be  normalized  to  unity  in  the  L2 
norm. 

For  any  body  the  lowest  eigenvalue  and  its  eigenfunction  are  always 

Ao  =  0  (2-15) 

V’o  =  Cl  (ci  is  a  normalizaton  constant)  (2-16) 

which  corresponds  to  uniform  current  injection  into  the  end  of  the  body  with 
no  current  sink;  such  a  current  has  A;o  =  0  and  therefore  is  constant  in  the 
^-direction;  it  rei^resents  common-mode  current  flow  through  the  body  to 
z  =  00. 


Since  we  are  interested  in  differential  mode  current  injection  we  will 
henceforth  assume  that  cq  =  0;  then  the  decay  of  the  spreading  current  in  2 
is  controlled  by  the  next  lowest  eigenvalue  Ai  >  0. 

ki  =  \/Ai,  (2-17) 

K  =  (8.69dB)A;i.  (2-18) 
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3  RESULTS  FOR  BODIES  OF  VARIOUS  SHAPES 


Example  1)  The  body  is  a  solid  rectangular  object,  so  the  cross- 
section  5"  is  a  solid  rectangle  of  dimensions  d  by  ic  in  the  x  and  y  directions 
respectively  (with  d  >  w).  Then 


—  Cl  cos(7ra:/d)  independent  of  ^ 
Ai  =  {n/df 
ki  =  Tx/d 
K  =  27dB/d. 


(3-1) 

(3-2) 

(3-3) 

(3-4) 


This  result  shows  a  remarkably  fast  fall-off  K.  For  instance,  to  achieve  an 
80  dB  fall-off,  we  need  only  go  in  the  long  direction  three  times  the  width. 

Example  2)  The  body  is  a  hollow  rectangular  box,  so  the  cross-section 
5"  is  a  hollow  square  of  dimensions  d  by  lo  in  the  x  and  y  directions  respec¬ 
tively  (with  d  >  w),  and  with  some  small  wall-thickness  t.  Then  to  a  good 
approximation 

V^i  =  Cl  cos(7rp/(d  -|-  w))  {p  is  distance  along  the  perimeter)  (3-5) 


=  {Tx/{d  +  w)Y 

ki  =  Tx/{d  +  w) 

K  —  27  dB /{d  +  tu) 


(3-6) 

(3-7) 

(3-8) 


independent  of  t  when  t  <C  ic. 

Example  3)  The  body  is  hollow  and  has  a  thin-walled  cross-section  of 
any  shape,  with  P  the  total  perimeter,  that  is,  the  length  of  C.  Then 


ki  =  2tx/P 
K  =  54dB/P, 


(3-9) 

(3-10) 
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approximately  independently  of  the  wall-thickness  t  P 


Example  4)  The  body  is  a  solid  circular  cylinder,  so  the  cross-section 
5  is  a  solid  disk  of  diameter  d.  Then  eigenfunctions  are  constructed  out  of 
Bessel  functions  J„(r)  and  angular  functions  exp(im0),  and 


If 


tjji  =  Cl  Ji(A;ir)  exp(i0)  (r  and  0  are  polar  coordinates)  (3-11) 


Ai  =  4{x'iq  +  l^)/d^  (xjo  is  the  first  root  of  J[) 
ki  =  4.19/d 
K  =  36dB/d. 


(3-12) 

(3-13) 

(3-14) 


Example  5)  The  body  is  an  “I-beam”,  so  the  cross-section  is  a  letter 
“I”  with  depth  d  (in  the  y-direction),  width  w  {<  d)  (in  the  x-direction),  and 
small  thickness  t.  Then  a  good  approximation  is 


t/’i  =  Cl  cos{ny/d)  independent  of  x 
Ai  =  {n/df 
ki  =  7r/d 

K  =  27dB/d. 


(3-15) 

(3-16) 

(3-17) 

(3-18) 


These  examples  are  summarized  in  Figure  2. 

Example  6  We  now  turn  to  bodies  constructed  of  discrete  components, 
such  as  lattices  of  metal  struts.  To  model  these  bodies  we  employ  infinite 
ladders  of  resistors,  as  shown  for  example  in  Figure  3(a).  Each  stage  consists 
of  one  transverse  resistor  (a  “rung” )  Ri ,  and  two  longitudinal  resistors  (the 
“stringers”)  i?2-  To  compute  the  properties  of  such  a  ladder,  use  a  well- 
known  trick  shown  in  Figure  3(b):  The  unknown  effective  resistance  Rx 
looking  into  the  end  the  ladder  can  be  determined  by  observing  that  adding 
one  more  stage  to  the  ladder,  leaves  it  invariant,  so 


1/ Rx  —  l/i?i  +  1/ {Rx  +  2i?2)) 


(3-19) 
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T  Solti 

J  ree+a*t^u,l«tr  l<[r  37 

i  C\ro<;S  5(^c^4iv^. 

If- VI  -A 


f 


l«-w 


*/«*») 

6+iAIV\  Wafi) 


K=  SV  ^/s/p 


H) 


Ccil'i  , 

Ctv'ole  Y\- '^(s  45 /d 


V^-vj  -M 


K=  s-p  ^/b/j 


Figure  2:  Cross  sections,  dimensions  and  values  of  K  for  the  examples. 
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Ct) 

Figure  3;  (a)  A  semi-infinite  resistor  ladder.  All  transverse  resistors  are  Ri. 
All  longitudinal  resistors  are  i?2-  (b)  The  effective  resistance  Rx  looking  into 
the  ladder,  which  can  be  determined  by  adding  one  stage,  which  gives  the 
same  Rx- 

which  is  readily  solved  for  Rx  as 

Rx  —  — Ri  T  R2  T  2i?i  R2  (3-20) 

For  instance,  if  all  resistors  are  equal,  Ri  =  R2,  the  fall-off  is  readily  shown 
to  be 


K  —  lldB/stage  ladder  of  equal  resistors  .  (3-21) 

To  achieve  higher  values  of  K  we  can  subdivide  the  /?2’s  and  add  more  rungs 
Ri  for  a  total  of  N  rungs  per  stage,  which  we  model  with  R^  =  R\/N.  We 
compute 
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N 

K 

1 

11  dB/stage 

2 

17dB/stage 

3 

21  dB/stage 

4 

24  dB/stage 

5 

27  dB/stage 

10 

39  dB  /  stage 

These  results  show  that  adding  more  “rungs”  increases  the  fall-off  substan¬ 
tially. 

Example  7  We  now  turn  to  3-dimensional  resistor  ladders,  as  in  Figure 
4(a).  Each  stage  consists  of  four  transverse  resistors  Ri,  and  four  longitudinal 
resistors  /?2-  How  shall  we  compute  this?  The  “ladder  trick”  does  generalize, 
but  becomes  cumbersome.  A  faster  computation  goes  as  follows.  The  ladder 
shows  a  four-fold  rotational  symmetry  when  we  look  along  it  from  its  end,  so 
the  eigenfunctions  of  current  distribution  must  be  proportional  to  exp(im0) 
for  m  =  (0, 1,  2)  where  0  is  angle  about  the  long  axis;  see  Figure  4(b,c,d).  As 
can  be  expected,  the  m  =  0  mode  is  common-mode  current  flow  along  the 
ladder  to  infinity,  and  can  be  ignored.  The  next-lowest  modes  are  the  two 
m  =  1  modes,  which  are  degenerate,  shown  in  Figure  4(c).  However,  in  these 
modes,  half  the  transverse  resistors  Ri  carry  no  current  because  the  voltage 
across  them  vanishes  by  symmetry.  Therefore  half  of  all  the  resistors  Ri  can 
be  chopped  out,  and  the  box  decomposes  into  two  2-dimensional  ladders  as 
in  Example  6).  Therefore  all  the  values  of  K  are  exactly  the  same  as  found 
above  for  Example  6). 

Again,  these  examples  apply  at  frequencies  so  low  that  inductive  effects 
can  be  neglected.  What  happens  at  higher  frequency? 
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Figure  4:  (a)  A  more  elaborate  semi-infinite  resistor  ladder,  with  a  box¬ 
like  stage.  All  transverse  resistors  are  Ri.  All  longitudinal  resistors  are  /?2. 
(b,c,d)  The  symmetries  of  the  normal  modes  of  this  ladder  are  like  exp{im4>) 
with  m  —  (0, 1,  2)  respectively. 
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4  FINITE-FREQUENCY  EFFECTS:  INDUC¬ 
TION 


As  frequency  a;  is  raised,  the  first  effect  to  become  important  is  induction 
of  electomotive  force  (EMF)  by  the  time-dependent  magnetic  field  B  (while 
charge  separation,  and  displacement  current  associated  with  time-dependent 
E  are  still  negligible).  Maxwell’s  equations  can  then  be  written  [2] 


V-J  =  0,  J  =  aE  (4-1) 

V  X  B  =  fiJ  (4-2) 

V  X  E  =  —dB/dt  (induction  equation)  (4-3) 

E  =  -V^-dA/dt  (4-4) 

B  =  V  X  A  (4-5) 

where  (4-6) 

A  =  vector  potential  with  V  •  A  =  0.  (4-7) 


Here,  all  fields  depend  on  time  like  exp{—iujt),  and  moreover  depend  on 
spatial  coordinates  {x,y,z);  fi  is  the  magnetic  permeability  of  the  material. 
Eliminating  E,  B  and  J  leads  to  a  Laplace  equation  for  4>,  and  a  diffusion 
equation  for  A 


=  0 

(4-8) 

V'^A  +  iixaijjA  =  — /r<7V$. 

(4-9) 

The  “skin  effect”  is  one  important  effect  governed  by  these  equations:  current 
tends  to  concentrate  near  the  surface  of  a  conducting  body,  in  a  layer  of 
thickness  roughly 

S  =  \j2/fiouj  “skin  depth”  (4-10) 
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with  fields  and  currents  falling  off  like  exp{—d/6)  with  distance  d  from  the 
surface.  When  this  is  the  main  effect  at  finite  frequency,  we  expect  Eq.  (3-10) 
above  to  still  apply  approximately,  with  the  “body”  replaced  by  a  thin  shell 
of  thickness  t  ^  6.  To  estimate  the  effect,  recall  that  skin  depth  is  of  order  a 
few  cm  in  common  metals  at  60Hz,  except  that  it  may  be  much  less  due  to 
permeability  of  ferromagnetic  materials.  Thus  for  frequencies  of  60Hz  and 
higher,  the  thin-shell  approximation  Eq.  (3-10)  applies. 


However  a  different  effect  may  be  much  more  important.  Unlike*  cur¬ 
rents,  magnetic  fields  can  traverse  free  space  outside  the  body;  different  parts 
of  the  body  can  therefore  be  coupled  by  mutual  inductance.  This  coupling 
can  be  expected  to  fall  off,  not  exponentially,  but  as  a  power-law  of  distance 
2.  (This  is  because  the  exciting  current  generally  has  some  magnetic  I’-pole 
moment,  and  the  magnetic  field  falls  off  like  (distance) from  it.  Tech¬ 
nically,  the  equation  for  A  has  an  unbounded  domain,  and  therefore  has  a 
continuous  spectrum  of  eigenvalues.) 


Let  us  estimate  the  inductive  effect  of  finite  frequency  /  in  a  rectangular 
body  of  cross-section  d  x  w  as  above.  If  the  current  source  has  a  magnetic 
dipole  moment  ~  Id^,  then  the  electric  field  induced  by  the  time-changing 
magnetic  field  propagating  through  free  space  is  roughly 

E^induct  ~  (4-11) 

where  uj  =  27r/  is  the  angular  frequency  and  is  the  constant  permeability 
of  free  space  (and  we  take  the  material  as  not  strongly  magnetic).  This  must 
be  compared  with  the  conductive  electric  field  inside  the  material 

E/conduct  ~  j  6^P(  kz')  .  (4-12) 

adw 


Their  ratio  is 

E^induct  dw  exp(kz)d? 

-  r\j  -  .  - 

E/conduct  ^ 


(4-13) 
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where  Eq.  (4-10)  has  been  used  to  eliminate  cu.  So  in  this  example,  inductive 
effects  dominate  spreading  current  at  modest  values  of  z,  when  S  is  no  larger 
than  the  transverse  size.  In  a  hollow  body  (or  “I-beam”)  of  thickness  t,  the 
corresponding  expression  is 

£^induct  dt  exp{kz)d'^ 

£^conduct  •  ^2  ^  -  ) 

and  inductive  effects  will  certainly  dominate  when  6  <  t. 

For  a  discrete  ladder  (as  in  Examples  6  and  7  above),  one  significant 
effect  is  that  each  component  now  has  a  combination  of  resistive  and  self- 
inductive  impedance.  When  this  is  the  main  effect  of  finite  frequency,  we  still 
find  exponential  decay  of  currents  along  the  ladder,  and  K  can  computed  in 
just  the  same  way.^ 

However,  once  again,  mutual  inductance  is  likely  to  dominate.  The 
mutual  inductance  of  distant  discrete  elements  generally  falls  off  as  a  power- 
law  in  stage  number  along  the  ladder,  not  as  an  exponential,  and  is  likely  to 
dominate  at  long  distances. 

For  both  continuous  and  discrete  bodies,  the  important  effect  at  high 
enough  frequency  is  therefore  likely  to  be  the  inductive  coupling  of  distant 
elements,  which  would  have  to  be  computed  by  different  methods  than  devel¬ 
oped  in  this  note.  At  yet  higher  frequencies,  charge  separation  and  capacitive 
coupling  will  eventually  become  important  too.  Finally,  at  RF  frequencies, 
displacement  current  becomes  important,  giving  radiative  effects. 


now  has  an  imaginary  part  which  describes  an  uninteresting  phase  shift  along  the 

ladder. 
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5  CONCLUSIONS 


We  have  shown  how  to  compute  the  fall-off  of  spreading  current  in  long 
conducting  objects,  for  differential-mode  DC  current  injection  at  one  end.  We 
have  also  computed  a  variety  of  examples.  The  results  show  that  the  fall-off  is 
remarkably  fast  in  many  examples.  Practical  applications  are  recommended 
when  the  control  of  spreading  currents  is  desirable.  AC  current  spreading 
will  differ  due  to  inductive  effects,  and  must  also  be  carefully  assessed. 
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